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Recently a covariant approach to cold matter universes in the zero-shear hypersurfaces (or longi¬ 
tudinal) gauge has been developed. This approach reveals the existence of an integrability condition, 
which does not appear in standard non-covariant treatments. A simple derivation and generaliza¬ 
tion of the integrability condition is given, based on showing that the quasi-Newtonian models are a 
sub-class of the linearized ‘silent’ models. The solution of the integrability condition implies a propa¬ 
gation equation for the acceleration. It is shown how the velocity and density perturbations are then 
obtained via this propagation equation. The density perturbations acquire a small relative-velocity 
correction on all scales, arising from the fully covariant general relativistic analysis. 


I. INTRODUCTION 

Density and velocity perturbations of cold matter universes are crucial to the understanding of structure formation 
in cosmology ||]. On scales well below the Hubble length, Newtonian theory may be used to analyze gravitational 
instability. However, current and upcoming observations and simulations are probing scales which are a significant 
fraction of the Hubble length, thus requiring a relativistic treatment [^. Bardeen’s gauge-invariant theory Q is 
most often used for relativistic perturbations. Various choices of gauge are possible in the theory. The zero-shear 
hypersufaces (or longitudinal) gauge sets up a frame of reference which emulates that of Eulerian observers in 
Newtonian theory, thus motivating the term ‘quasi-Newtonian’. Since the frame is non-comoving, there are relative- 
velocity effects on the density perturbations, and subtle issues arise in dealing with these effects, as pointed out 
recently by van Elst and Ellis |^. In order to resolve these problems in a fully gauge-invariant way, a covariant 
approach may be adopted. 

The covariant and gauge-invariant approach to perturbations was developed by Ellis and Bruni on the basis 
of Hawking’s paper . It has various advantages over the non-covariant gauge-invariant theory (see for further 
discussion). One advantage is that all quantities have a direct and immediate physical or geometric meaning, and 
no nonlocal decomposition into scalar, vector or tensor modes is required. A second advantage is that the covariant 
approach provides a natural and transparent setting to search for integrability conditions which may arise from 
constraints. These general relativistic constraints and the consequences that follow from their evolution, are often not 
made explicit. As a result, crucial general relativistic effects can sometimes be obscured or missed. 

Covariant consistency analysis of constraints was developed by Maartens |^, building on the methods of Lesame 
et al. ]l^ , in a form applicable to both the nonlinear exact theory and the case of linearized perturbations. (Further 
developments of the approach are given in 011 .) Applications of a covariant approach to ‘silent’ universes [00, 
to nonlinear gravitational radiation [^7|j^] , and to non-accelerating fluid models reveal the existence of crucial 

integrability conditions. The covariant characterization of scalar, vector and tensor perturbations also relies on such 
an approach [0-0. 

Van Elst and Ellis [^ use a covariant consistency analysis first in nonlinear quasi-Newtonian cosmologies and then in 
the case of covariant linearization about a Friedmann-Lemaitre-Robertson-Walker (FLRW) background. They show 
that the nonlinear models are likely in general to be inconsistent, except for special cases such as FLRW solutions. 
Inconsistency of the nonlinear models is not surprising, since one is demanding for all possible dynamical evolutions 
of matter that there exists a shearfree and irrotational congruence, which in particular forces the magnetic part of 
the Weyl curvature to vanis h [0|. This rules out gravitational radiation |^, and thus leads to severe restrictions 
on the gravitational held [ p^27]j28( | . In the linearized case, one might expect that all the problematic terms which 
arise in evolving the constraints would be removed. It is implicitly or effectively assumed in standard non-covariant 
perturbation theory that there are no integrability conditions arising from the zero-shear hypersurfaces gauge. 

However, it turns out that the linearized models are not in general consistent. Van Elst and Ellis introduce an 
ansatz for the evolution of the gravitational potential, which they motivate by a discussion of the lapse in ADM- 
type approaches. Using this ansatz, they hnd an integrability condition in the linearized models. The integrability 
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condition is satisfied by a particular value of the constant parameter in the ansatzj^ The reason that the integrability 
condition is not revealed in some non-covariant treatments is probably either an implicit assumption that only evolution 
equations need be considered, or a gauge-dependent approximation that effectively neglects the velocity and removes 
the constraints (see [§ for further discussion of this point). 

In this paper, extensions of the results of on linearized quasi-Newtonian models are given. The main result is 
the determination of the velocity and density perturbations. A crucial part of the analysis is the combined use of the 
comoving (‘Lagrangian’) and quasi-Newtonian (‘Eulerian’) frames. 

Section II summarizes the necessary covariant equations and methods, with some details given in appendices. 
Section III uses a transformation to the comoving frame in order to show that the quasi-Newtonian models are a sub¬ 
class of the linearized silent models. This is the basis for a simple and direct approach to deriving the integrability 
condition in general, i.e. without introducing any ansatz for the gravitational potential. The general integrability 
condition reduces to the special form given in when one imposes their ansatz for the gravitational potential. Their 
ansatz is also generalized. Furthermore, another integrability condition is derived by considering spatial consistency of 
the constraints. The van Elst-Ellis solution of the first integrability condition is shown to reduce the second condition 
to an identity. 

The main result however, follows from the fact that the van Elst-Ellis solution itself implies a crucial propagation 
equation for the 4-acceleration. This propagation equation then leads to an equation determining the velocity pertur¬ 
bations. The equation is scale-independent, so that velocity perturbations have an effect on all scales. The velocity 
perturbation equation is readily solved analytically for a flat background. 

In section IV, the density perturbations are found by a direct and simple approach. The complicated calculations in 
1^ of the energy flux (momentum density) source term in the density perturbation equation are by-passed by consid¬ 
ering the density perturbations in the comoving frame. Furthermore, the equation is solved for a flat background. The 
covariant correction to density perturbations that arises from relative-velocity effects is a simple comoving divergence 
term, which can be found from the velocity perturbations. This correction affects all scales, although it is rapidly 
dominated by the usual solution. The correction to the growing mode is consistent with the result of Takada and 
Futamase [Q, who use non-covariant theory, but not in the zero-shear hypersurfaces gauge. 

Concluding remarks are made in section V. 


II. COVARIANT EQUATIONS 


Given a choice of 4-velocity field (with = — 1), the Ehlers-Ellis approach employs only fully covariant 

quantities and equations with transparent physical and geometric meaning. The quantities are split into spacetime 
scalars and spatially projected tensors, while the equations split into evolution equations along and constraint 
equations involving only spatial covariant derivatives. These equations arise from the Ricci identity for and the 
Bianchi identities, with Einstein’s field equations incorporated via algebraic replacement of the Einstein tensor by the 
energy-momentum tensor Tab- The covariant linearization of the equations is the basis for the Ellis-Bruni perturbation 
theory [^. Integrability conditions, at both the nonlinear and linear levels, arise from investigating the derivatives 
of constraint equations (including additional conditions that may be imposed by physical or geometric assumptions), 
using covariant differential identities and the evolution equations. The streamlined and developed version of the 
Ehlers-Ellis formalism given by Maartens [|| (see also |^,^,|2|) greatly facilitates such investigations, especially 
by making explicit the irreducible quantities and derivatives, which significantly simplifies the equations, and by 
developing the covariant identities which these quantities and derivatives obey. 

The projection tensor hab = gab + UaUb, where gab is the spacetime metric, and the projected alternating tensor 
£abc = Vabcd.u'^, where r]abcd = — [a^i^b^"^d] IS the spacetime alternating tensor, are the basis for covariant 

irreducible splitting of tensors and derivatives. Projected rank-2 tensors Sab are split into a scalar trace, a projected 
vector spatially dual to the skew part, and a projected symmetric tracefree part: 


Sab = \Sadh‘^‘^hab+£abcS‘^ + S. 


(ab) 5 


where Sa = ^SabcS^’”^^ = S(^a} = habS^ and iS'(ah) = [h(a'^hb)'^ — ^h'^‘^hab]Scd- Covariant time and spatial derivatives are 
defined by 

S‘^ 'b- = u'^VcS^ 'b- , DcS'“ '6- = hjh'^d ■ • ■ /lb® • ■ • VfS'^ 'e... , 


^ For other values of the parameter, only very special models appear to be consistent . 


2 



and then the covariant spatial divergence and curl are 

divF = D“V;, {divS)a=B^Sab, (1) 

curl K = SabcOV^ , curl S'ab = ecd(aD°S'b)‘^. (2) 

The dynamic quantities are the energy density p, the pressure p, the energy flux qa = q{a) i and the anisotropic stress 
TTab = 7r(ob), so that 

Tab — P'^a'^b “t" P^ab “t“ ‘^q(^a^b) '^ab ■ 


The kinematic quantities are given by 

^b^a — Db^a -^a^b ; Db^a — g^^ab “f ^ab “t“ £abc^ : 

where 0 = D“ua is the expansion, Aa = iia = ^(a) is the 4-acceleration, uja = —^curlrta = W(a) is the vorticity, and 
CTab = D{a'*^b) is the shear. Finally, the gravito-electromagnetic fields are 

Tab — CacbdU U — ^[ab) ; ^ab — 2^acdC be^ — ^[ab) i 

where Cabcd is the Weyl tensor, which represents the locally free gravitational field The FLRW background 

is then covariantly and gauge-invariantly characterized by: dynamics - DaP = 0 = DaP, qa = 0, iTab = 0; kinematics 
- Da0 = 0, Aa = 0 = Wa, cdab = 0; gravito-electromagnetic field - Eab = 0 = Hab- 

In this paper only the linearized quasi-Newtonian cosmologies are considered, since the main focus here is on 
perturbations and structure formation. The covariant linearized evolution equations in the general case are 


P =-(P + P)0 - divg, (3) 

0 = -^0^ - 5 (p -I- 3p) -I- div A , (4) 

Qa = -dHqa - (p + p)Aa - DaP - (div 7r)a , (5) 

OJa = -2Huja - icurl Aa , (6) 

dab — (7 ab Eab T 2 '^ab 4“ ^ {a-^b) : (^) 

Eab = -3HEab + cm\ Hab - ^dab “ ^(p +p)CTab “ ^D(agb) “ \HTTab , (8) 

Tab = -3TTab - curl iJab + ^CUrlTTab , (9) 

and the linearized constraint equations are 

= divw = 0 , (10) 

C^a = (divcr)a -curlwa - fDa0-t-ga = 0, (11) 

0 ab = curl (Tab 4“ D^a^b) Eab — 0 , (4^) 

C^a = (divF;)a 4- d(div7r)a - ^I)aP + Hqa = 0, (13) 

C^a = (div H)a + dcurl qa - {p + p)uJa = 0 , (14) 

where H = a/a is the background Hubble rate, related to the background values of p and p by 

P = 3(^T2 + ^^ , T + T" = -i(p + 3p), (15) 


with /C = 0, 4il the curvature index. The differential identities that are needed for investigating consistency and 
deriving perturbation equations are collected in Appendix A. 

If another 4-velocity is chosen, the corresponding kinematic, dynamic and gravito-electromagnetic quantities 
undergo transformations. For completeness, and since they do not appear elsewhere, the exact nonlinear form of 
these transformations is given in Appendix B. Only the linearized form of the egressions is required below. 

The covariant characterization of quasi-Newtonian cosmologies is as follows g: they are almost-FLRW dust uni¬ 
verses, with a congruence of observers whose 4-velocity field u°‘ is irrotational, shearfree a nd nonrelativistic relative 
to comoving observers. The comoving 4-velocity is given by the linearized form of Eq. (Bl): 




( 16 ) 
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where is the nonrelativistic relative velocity, which vanishes in the background. The models are thus defined in 
the quasi-Newtonian frame by 


(dynamics) p = 0, Qa = pva , TTab = 0 , (17) 

(kinematics) oJa = 0 , aab = 0 , (18) 


as shown by the linearized form of equations dUli-dl). Thus quasi-Newtonian cosmologies are irrotational, shearfree 
dust spacetimes with energy flux (momentum density) arising purely from a particle flux in the quasi-Newtonian frame 
relative to the comoving frame. Note that the isotropic and anisotropic stresses that arise from relative motion are 
second order in Ua, as given by equations (p^) and (BIO). 


The gravito-magnetic constraint equation (1^), together with Eq. 


shows that 


Hab = 0. (19) 

Thus there is no gravitational radiatio n ||l5| , which further justifies the term ‘quasi-Newtonian’. In addition, the div 
H constraint (0), together with Eq. (|18|), shows that qa is irrotational, and thus so is Va- 

curl Uq = 0 = curl ( 7 q . (20) 


Since the vorticity vanishes, it follows (see Q) that Va 


Daip, where the velocity potential "0 is determined below. 


III. INTEGRABILITY CONDITIONS 

In irrotational dust models with vanishing energy flux and anisotropic stress, the constraint equations = 0 evolve 
consistently with the evolution equations, even at the nonlinear level, in the sense that (see also UHl) 

= F^b + G^Ba , 

where F and G depend only on the kinematic, dynamic and gravito-electromagnetic quantities (and not their deriva¬ 
tives). If one imposes the ‘silent’ constraint ([l^, then the nonlinear models are generically inconsistent, but the 
linearized models are consistent p^ , p^ . A very simple approach to the integrability conditions for quasi-Newtonian 
cosmologies follows from showing that these models are in fact a sub-class of the linearized silent models. This can 
be seen by transforming to the comoving frame. 

Linearizing the expressions in appendix B for the case where and iF are any frames in nonrelativistic relative 
motion, one finds for the kinematic quantities 

0 = 0-1- div V , 

Aa — Afi -\- Vd -t- tiVa , 

COa=LOa- |curlWa , 
dafe — ^ab T 19(a^b) 5 

for the dynamic quantities 


p = P, P = P, qa=qa- {p+ p)Va , ^ab = T^ab , (25) 

and for the gravito-electromagnetic field 

Eab = Eab , Hab = Hab ■ (26) 

Eor the quasi-Newtonian frame and iF the comoving frame, it follows from the equations (0)-@ that 

p = 0, go = 0, TTob = 0, (27) 

Ao — 0 , U)a — 0 t O'ob — I-l(a^6) ; (28) 

Eab = Eab , iiab = 0 . (29) 


Equations (|^)-(|^) constitute a covariant characterization of linearized silent universes, except that the shear takes a 
special form. Thus quasi-Newtonian models are linearized silent models with a special form of shear, and integrability 
conditions arise only from the restriction on the shear. 


( 21 ) 

( 22 ) 

(23) 

(24) 
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It is now convenient to return to the quasi-Newtonian frame, where the restriction on the shear is that it vanishes. 
Integrability conditions arise directly from the fact that the shear propagation equation is turned into a constraint, 
i.e. Eab = D(o^b)- This can be simplified, using curlA^ = 0, which follows from the vorticity propagation equation 
(^, and identity (^). Thus 


Aa = DaV?, (30) 

where ip is the covariant relativistic generalisation of the Newtonian potential. Then the shear constraint becomes 


^ab — Efib — 3 • 


(31) 


In summary, the only independent new constraint is (31), and any conditions that follow from its derivatives. What 
is happening here is that the consistent evolution of the basic constraints {A = 1, 2, • • •, 5) is not affected by 
introducing a new constraint. It is the new constraint E which leads to integrability conditions. The freedom in the 
gravito-electric field is clearly central to the consistency of the silent models, and conversely, it is the longitudinal 
condition (^^ on that field which produces integrability conditions in the quasi-Newtonian subcase. 


A. Time evolution 


The time derivative of Eq. (M) follows from the gravito-electric propagation equation (H) and the identity 


{D{aD6)‘/^} = D(aDb>V3+ - 2iI)D(aDf,)¥J, 

which is proved using identities (^) and (^) and Eq. (p0|). It follows that 

Eab = —^HEab ~ — [if + H) — D^aDfe) {ip + §0) I (32) 

on using Eq. ©■ Thus Eab evolves consistently if 

D{aD;,) {ip + |0) -I- -I- |0) = 0 . (33) 

This is the first integrability condition in quasi-Newtonian cosmologies. It represents an extension of the condition 
derived in |^, since there a particular ansatz is assumed a priori for ip, i.e. 

ip = Q!0 , 

where a is a constant parameter. Choosing a = i.e. 

(^+|0 = O, (34) 

it is clear that the integrability condition d^ ) is reduced to an identity. 

The equation ( ^4|) will be adopted here, since it is covariant, has a clear geometric motivation (as given in j^), and 
guarantees consistent evolution of the gravito-electric constraint. Before proceeding with the van Elst-Ellis solution, 
it is interesting to ask whether it may be generalized. The integrability condition ( |33[) may be rewritten as 

D<„Df,){e‘"((p+i0)} = O. (35) 

This shows clearly how the van Elst-Ellis solution may be generalized to 

-I-10 =/3e“‘^ , (36) 

where Da/3 = 0, i.e. /3 is an arbitrary background scalar. There does not appear to be any advantage in adopting the 
generalized solution (^). It is not clear whether more general solutions of the condition may be found. 

What are the immediate consequences of the van Elst-Ellis solution to the integrability condition? Firstly, the time 
evolution of (|3^) itself leads to the covariant modified Poisson equation 

I)‘^ip= ^p-{Sip + Sif) , ( 37 ) 

after using the Raychaudhuri equation (|^). This equation governs the relativistic gravitational potential for a given 
energy density. 
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Secondly, one can get a crucial evolution equation for the 4-acceleration |Q. Such an evolution equation is not 
present in the set of general evolution equations. It aris es via the shearfree condition, as a consequence of Eq. 
Taking the gradient of Eq. (^), and using identity (AS) and the div a constraint d^), one gets! 


Aa + 2HAa — —\pVa 


(38) 


Now there is also an evolution equation for Va 


Va + HVa = -Aa , (39) 

as follows from the conservation equations (^) and (||), or from the comoving frame equations ( p^ and (|^). This is 
just the relativistic generalization of the Newtonian equation for relative acceleration: 


dv 

dt 


—V(p. 


The coupled evolution equations (|^ ) and (|^ ) may be decoupled to produce second order equations in either quantity. 
For Va- 


Va + 3HVa 



Va = 0, 


on using Eq. (p^), while for Aa'. 

Aa + (>HAa ^ - Aa = 0 . 


(40) 


(41) 


These equations may be solved to find the velocity perturbations Va and the 4-acceleration Aa- Since there are no 
spatial derivatives in Eq. (p^, the velocity perturbations are independent of scale. 

For a flat background, with AS = 0 and H oc equation (p0|) is readily solved analytically: 


= Ai+) + Ai-) 




(42) 


where (-I-) and (—) denote the growing and decaying modes, and = 0. Using Eq. (|^, the solution for the 
acceleration follows as 


A„ = -A(+)a-i + |Ai-)a-^/^ (43) 

Finally, one can derive another equation for ip using Eq. (^), which implies, together with Eq. (pT|), that 

2 

Va =- ^a‘P ■ (44) 

P 

Then using Eq. (^) in Eq. (po|), it follows that 

D„^ + 3iID,(^-i(p-3El2)D„(^ = 0. (45) 

Note that the background energy conservation equation and Eq. (^) determine the velocity potential in terms of the 
gravitational potential: 

/ 2 \ 3 

Va = Da'0 where ip = — { -r a p. 

\Po%J 


Note that the same evolution equation (B® is obtained if the generalized solution (pm is used, with [ 3^0 — Da/3. 
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B. Spatial consistency 


Spatial derivatives of the new constraint £ determine whether it is consistent with on an initial spatial surface. 
Taking the curl of Eq. (^) produces no restrictions, since cnrlEab = 0 by the gravito-magnetic propagation equation 
a nd c onstraint equation (p^, and since curlD^^jD^^c/j vanishes identically. This follows from using identities ( [^ ) 
and (A6): 


curlD^aDf,)(/? = curlD(aDt,)(/3 = curlDaD^v? 


The divergence however gives a nontrivial condition. First, one needs an identity for the divergence of the distortion: 


+ iD„ (div A) + i (p - , 


(46) 


which holds for any projected vector Aa, and follows from identity ([A6| ). 
that 


Using equations (|30|) and (|A2|), it follows 


D^D(,D6)(p = |D, (dV) + I (p - 3i?') . (47) 

Now use equations and (^ to get 

(div£)^ = C\ - HC\ + iD,p - |i7D,0 - |D, (dV) - § (p - 3i7") D„(p. 

Thus there is a second integrability condition arising from Eq. (^), i.e. 

D„p - 2i7D,0 - 2D„ (dV) - 2 (p - D,(p = 0 . (48) 

In general, this appears to be independent of the first integrability condition (|^). However, if one uses the van 
Elst-Ellis solution (|34[ ) of the first condition, then the second condition is identically satisfied. This can be seen as 
follows. Taking the gradient of Eq. (^) and using Eq. (^, one finds that the second integrability condition ( |4^ ) 
becomes 


2 [Ba(p + 3i7D,(p - |(p - 3772) ^ q ^ 

which is identically satisfied by virtue of Eq. (^) . 


(49) 


IV. DENSITY PERTURBATIONS 

In 1^ the density perturbation equation is found in the quasi-Newtonian frame, which requires incorporating a 
complicated energy flux source term. A simple alternative is to work in the comoving frame, which leads directly to 
solutions of the density perturbation equation in terms of the velocity perturbations. 

The covariant density perturbation scalar in the comoving frame is 

5 = aD 5a where 5a = —— ■ (50) 

P 

Using Eq. (|^ ) and the identity 

Da/ = Da/ + fVa , 

it follows that to linear order the density perturbation scalars in the comoving and quasi-Newtonian frames are related 

by 

5 = I -f 3o2i7D“ua . (51) 

Covariant time derivatives in the comoving and quasi-Newtonian frames agree to linear order: lEVaS = S. Thus 
5 satisfies the standard equation (with 1C = 0) for dust Q 
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(6)" +2H{Sy = 0. (52) 

Since the solution 6 of Eq. (^) is well known, one can use Eq. ( pl| ) to write down the density perturbations in the 
quasi-Newtonian frame: 

(5 = A(+) aA(-) , (53) 

where = 0 , and 7 = is a constant, with oi = and 


V = oD’^Va 


is the comoving divergence of the peculiar velocity held. This quantity encodes the scalar contribution of velocity 
perturbations to density perturbations. Using the solution ( ^ for Uo, one hnds that 

V = r(+) 0^/2 p(-) q- 2 p(±) ^ . (54) 


It follows from the identity (A4) that E^ = 0. Then using Eq. (p4|) in Eq. (p^ gives 


5 = A(+) a -1- A(-) a-3/2 -h 7 r(+) -h 7 r(-) 


(55) 


The correction to the standard solution affects all scales. The growing mode is shifted by an amount that is a comoving 
constant, which will be dominated by the term that grows as a. A similar feature arises in a non-covariant Lagrangian 
perturbation analysis (see [^). The decaying mode (not considered in |^) is corrected by an amount which dies away 
more rapidly, and so becomes negligible. 


V. CONCLUDING REMARKS 


The covariant approach to quasi-Newtonian models developed by van Elst and Ellis ^ has been applied and 
extended here, in order to derive and solve the equations governing density and velocity perturbations. The quasi- 
Newtonian zero-shear hypersurfaces gauge (or longitudinal gauge) implicitly involves integrability conditions - i.e. 
it incorporates a dynamical condition, and is not pure gauge. This has been overlooked in some non-covariant 
treatments, effectively amounting to gauge-dependent assumptions about relative-velocity effects. A fully covariant 
general relativistic treatment has uncovered these integrability conditions and their implication for the perturbations. 
The perturbations have been determined, making explicit the relative-velocity effects. 

The integrability condition found in 1^ was generalized, via a simple approach that showed how the quasi-Newtonian 
models are in fact a sub-class of the linearized silent models. The generalized first integrability condition (|^) natur ally 
leads to the generalization (|^) of the van Elst-Ellis solution (Q). Furthermore, a second integrability condition ( |4^ ) 
was found by investigating spatial consistency. The integrability conditions were shown to be crucial for determining 
the perturbations, since the solution ( ^^ of the integrability conditions implies a propagation equation ( |^ ) for the 
4-acceleration. In turn, this propagation equation leads to the velocity perturbation equation ( [l0[), and the solution 
( ^^ was given for a flat background. By transformii^to the comoving frame, a simple relation (^iD was derived for 
the density perturbations, which were given by Eq. (|53| ) in terms of the velocity perturbations and the standard dust 
solution for a flat background. The density perturbations were given analytically by Eq. (55), which showed how 
relative-velocity effects produce small corrections to both the growing and decaying modes on all scales. 

These results underline the importance of general relativistic constraint equations and the integrability conditions 
that can arise from imposing various physical and geometric assumptions in cosmology. A covariant approach [^0|J^] 
avoids many of the intricate problems arising from gauge-dependent approaches. The improved covariant formalism 
of [^, summarized in section II and supplemented by the identities of appendix A and the new results in appendix 
B on transformations of the covariant quantities, has been central to deriving these results. Similar methods can in 
principle be applied to investigate the consistency and underlying implications of other gauge choices or of various 
approximations (such as the Zel’dovich approximation) in general relativistic cosmological perturbations. This is a 
subject of further research. 





APPENDIX A: LINEARIZED IDENTITIES 


Useful differential identities |Q: 


curlBaf =-2fuja , (Al) 

(D„/) = D„ (DV) + Hp- 3^") Da/ + 2/curl , (A2) 

(D„/)- = bJ - HBJ + fAa , (A3) 

{BaSb...y = BaSb... - HBaSb ..., (A4) 

(D V) ■ = D V - 2HD^f + /div A, (A5) 

D[aD,] 14 = i (3iL^ - p) V[a/l6]c , (A6) 

= I {3H^ - p) S[,^-hby ^^, (A7) 

divcurlU = 0 (A8) 

(div curl S)a = ^curl (div S)a , (A9) 

curl curl Vo = Da(divU) - + | (p - 3-ff^) Va , (AlO) 

curlcurlS'ab = |D(„(divS'),,) - D^Sab + (p - 3iS^) Sab , (All) 


where the vectors and tensors vanish in the background, Sab = S/^ab), and all identities except & are linearized. 
(Nonlinear identities are given in [^[p^J^.) 


APPENDIX B: TRANSFORMATIONS UNDER CHANGE OF FRAME 

Change in 4 -velocity: 

Ua=l{ua + Va) where 'CqU® = 0, 7 = (1 -'(;^)“D 2 . (Bl) 

The following algebraic relations are needed: 

9 ab — ^ab — ^ab ^a^b 7 

hab — T 'y UaUb T ^^[a^b) T 7 

‘^abcd — ^'^[a^b]cd ^^ab[c'^d] ~ ^'^[a^b]cd ‘^^ab[c'^d] 7 

^abc — 'y^abc T T T Uc^abd^ ^ 7 

Cab‘=^ = 4 {u[aU^’^ + /l[a''=} Eb]‘^^ + 2 eabeU^^ 

= 4 + 2u^aHb]ee^'^^ , 

together with the decomposition 

^bVa = -Ub {w{a> + AcV'^Ua] + Ua {56'^;, + CTbcV^" + [u!, v]b} 

+ i (divy) hab - ^EabcCurlu^ + D^^yb) , (B 2 ) 

where [W, V]a = £abckU^U‘^. Then the following exact nonlinear transformations may be derived. 

Kinematic quantities (using Va 7 = aVb)'- 

0 = 7© + 7 (div V + A°'Va) + 7^IU , (B 3 ) 

Aa = J^Aa + 7^ {■y(a) + §OVa + dabV^ - [uj, -uja + + UbcV^v'') Ua 

+ ^{divv)va + i[r),curlz)]a + ZJ^’D^b^^a)} + l^W{Ua + Va) , (B4) 

Cba=l^ {{I- \v^) IVa - icurlz;^ + ^Vb (2a;^ — curlu^) Ua + ^VbUj^Va 

+ ^[A,v]a + ^[v,v]a + ^SabcCr''dv‘'v‘^} , (B 5 ) 

dab = JCTab + 7(1 + 'l'^)U(a<^b)cV‘' + 7^A(a [z)b) + v'^Ub)] 
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+ 7D(o1'&) - \hab [Acv'" + 7^ (VF - Vcv'')] 

+ 'y^UaUb [tTcdf^'y'* + ^v‘^Acv‘' - V^v’^Bi^^Vd) + (7"^ - - 1) w] 

+ l^U{aVb) [AcV‘' + dcdV^'v'^ - VcV‘' + 2 ^^ (7^ - i) W] 

+ ll^VaVb [divi; - A^v^^ + 7^ (87^ -l)W]+ 'y^v^aib) + v^l^U(aV{b)) 

+ 7^'y(aCr&)cl^° - {l^b) + V^Ub)} + 2^^^ v'^D l^c.V (a) {vb) + M&) } , 


where W = + |z;^divz; + 

Dynamic quantities (compare [Q): 


Vd)- 


p = p + ^‘^ [v'^{p+p) - 2qaV’^ + ■KabV'^v'’] , 

V=V+\l'^ [u^(p +p)- 2qaV°- + TTabV°‘v'^] , 

go = iqa - ITTabV^ - 7 ^ [{p + p) - 2qbv" + TTbcV^Va 

- 7 ^ [v^{p+p) - (1 + V^)qbV^ + TTbcU^U®] Ua , 

^ab = TTob + 27^?;‘'7re(a {wb) + Vb)} “ 2w^7^g(aWb) “ ‘^I'^QiaVb) 
-^ 7 ^ [v'^ip + p) - 2 gcU^ + 7rcdv''v‘‘] hab 

+ \'y'^ [2u^(p + p) - 4u^gcf^ + (3 - v‘^)T:cdv''v'^] UaUb 

+ § 7 ^ [2v‘^{p + p) - (1 + 3v'^)qcV^ + 27rcd'uV] Uf^aVb) 
+ [(3 - v‘^)ip + p) - ^cV'' + 2TTcdV''v‘^] VaVb ■ 


Gravito-electromagnetic field: 

Eab = 7^ {(1 + v"^)Eab + v‘^ [2ej,d(a^b)^ + “^E^f^g^Ub) 

-\- {UgVb “t“ hgb^EgdV ‘^Eg^gVb) ‘2ll^g£b)cdE V j" , 

Hab = 7 ^ {(1 + V^)Hab + v‘^ [—‘^£cd{aEb)'^ + 2 ^fc(a'*^b) 

“t“ {VaUb hab^HgdV 2i7c(a^b) ‘^V(^a^b)cdE VaJ } ■ 


(B6) 

(B7) 

(B8) 

(B9) 


(BIO) 

(Bll) 

(B12) 
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